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1.  INTRODUCTION 


Cross-lagged  panel  studies  are  statistical  studies  in  which 
two  or  more  variable  are  measured  for  a  large  number  of  subjects 
at  each  of  several  waves  or  points  in  time.  The  variables  divide 
naturally  into  two  sets  and  the  primary  purpose  of  analysis  is  to 
estimate  and  test  the  strength  of  the  relationship  between  the 
sets.  Such  studies  are  found  in  the  mainstream  of  social,  behav¬ 
ioral  and  business  research  (e.g.  Crano,  Kenny  and  Campbell 
(1972),  Greenberg,  Kessler  and  Logan  (1979),  Eaton  (1978), 
Polachek  and  McCutheon  (1983)  and  Frey  (1984)) .  This  paper  con¬ 
tributes  to  these  studies  by  developing  and  applying  procedures 
for  detecting  the  presence  of  serial  correlation  in  the  error 
structure  of  the  regression  models  used  in  such  studies. 

Methods  of  analysis  for  cross-lagged  panel  studies  have 
evolved  over  the  last  twenty  years  (Kessler  and  Greenberg  (1981) 
provide  an  excellent  review  of  this  development).  Early  methods 
used  correlation  statistics  and  were  motivated  by  the  work  of 
Donald  Campbell  (1963)  ,  which,  in  turn,  was  motivated,  by  the 
seminal  work  of  Paul  Lazarsfeld  (1948)  on  the  analysis  of  panel 
studies  involving  discrete  variables.  The  most  significant 
advancement  in  the  methodology  of  cross-lagged  panel  analysis  has 
been  the  development  of  a  regression  approach.  This  approach 
treats  the  cross-effects  as  parameters  in  one  or  more  regression 
models  and  then  applies  standard  regression  procedures  to  estimate 
and  test  these  parameters  (for  example,  see  Peltz  and  Andrews 
(1964),  Duncan  (1969),  Heise  (1970),  Hannan  and  Young  (1977),  and 
Rogosa  (1980)).  Although  the  cross-lagged  panel  study  includes 
observations  over  time,  for  the  most  part,  this  approach  has 
assumed  that  the  errors  of  the  model  are  independent  across  waves. 

Mayer  (1984a;  1984b)  extends the  regression  approach,  by 
incorporating  the  cross-effects  as  parameters  in  a  multivariate 
regression  model  and  develops  procedures  to  estimate  and  test 
these  parameters.  He  considered  both  the  model  with  independent 
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errors  end  Che  model  with  serially  correlated  errors.  Xn  this 

/ 

paper  we  extend  the  applicability  of  this  results  by  considering 
the  problem  of  detecting  if  the  multivariate  regression  model 
requires  the  serially  correlated  error  structure  as  opposed  to  the 
independent  error  structure. 


2.  PRELIMINARY  NOTIONS 

To  formalize  the  multivariate  regression  structure  of  the 
cross-lagged  panel  model  let  z  ■  (x^;  ^c) '  be  the  measure¬ 
ments  made  on  the  k  variables  for  the  ith  replication  at  the  tth 
wave  (i  -  1,  ....  n;  t  »  0,  ....  T)  where  and  are  sub¬ 
vectors  of  p  and  q  measurements,  respectively.  Let  Z  ■  (X  ;  7  ) 
be  the  n  x  k  response  matrix  that  has  z  as  the  ith  row.  The 
regression  structure  is 


(2.1) 


It  *  fcc-1  l*  It  C  "  1 . T 

where  B  is  a  k  x  k  matrix  of  unknown  regression  coefficients  that 
does  not  depend  on  t;  and  Ec  is  an  unobserved  random  error  matrix 
with  structure 


It  ■  St-i  i  -  It  =-1 . 1 


(2.2) 


where  £  is  a  k  x  k  unknown  matrix  of  autoregressive  parameters  and 
£c  is  an  unobserved  random  matrix  with  rows  which  are  independent 


Gaussian  vectors  with  common  mean  0  and  common  covariance  matrix 


A  that  does  not  depend  on  t. 


The  rows  of  £c  are  thus  independent  Gaussian  vectors  with 


common  mean  0  and  common  covariance  matrix 


£c  “  Vtc-l  £  *  k  C"1 . T> 

If  we  let  0  be  the  covariance  matrix  of  Z  . 

ch«n 

^t  "  ®  tt-i  ®  +  ®'^t-i  ®  *  k  ^t-i  i  +  k 


(2.3) 


(2.4) 
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To  complete  the  specification,  we  describe  the  behavior  of 
the  response  matrix  and  error  matrix  at  the  initial  wave.  Since 
the  initial  observations  are  generated  by  the  same  sampling  scheme 
as  the  other  observations  we  assume  the  initial  response  matrix  ^ 
is  random  with  rows  that  are  independent  Gaussian  vectors  with 
common  mean  £  and  common  covariance  matrix  ^ .  Likewise,  we  let 
the  rows  of  the  initial  error  matrix  ^  be  Independent  Gaussian 
vectors  with  common  mean  ^  and  common  covariance  matrix 

The  primary  goal  in  the  research  approach  is  to  estimate  the 
regression  parameter  matrix  £  and  to  test  the  hypothesis  of  no 
effects  (£  ■  £  ■  0)  and  the  hypothesis  that  the  responses  in  ^ 
do  not  affect  the  responses  in  7c>  and  conversely,  the  responses 
in  Yt_^  do  not  affect  the  responses  in  -  the  hypothesis  that 
£  and  £  are  block  diagonal.  The  method  of  estimation  and  testing 
depends  heavily  on  whether  serial  correlation  is  present,  l.e., 
whether  ^  For  example,  if  serial  correlation  is  present, 

but  Ignored,  the  estimator  of  £  derived  by  applying  maximum 
likelihood  (to  the  wrong  model)  is  not  consistent  in  the  number 
of  replications  (see  Section  4). 

Before  presenting  our  tests  we  comment  on  the  appropriateness 
of  existing  procedures  for  detecting  the  presence  of  serial  cor¬ 
relation  in  regression  models  involving  observations  over  time. 

The  most  widely  used  procedure,  the  Durb in-Watson  procedure,  is 
not  appropriate  for  the  cross-lagged  panel  study  because  it. is  not 
designed  for  use  with  models  that  contain  only  lagged  endogenous 
predictors  and  is  designed  for  use  with  univariate  regression 
models  (Malinvaud  (1980),  Johnson  (1972)).  The  most  notable 
extensions  of  this  procedure  are  appropriate  for  models  with 
lagged  endogenous  predictors  but,  again,  are  not  designed  for 
multivariate  regression  models  (Durbin  (1970),  Judge,  et.  al. 
(1980)).  Finally,  we  know  of  one  extension  of  these  procedures  to 
a  procedure  chat  is  appropriate  for  multivariate  models  (Guilkey 
(1975));  but  it,  like  all  of  these  procedures,  guarantees  chat  the 
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advertised  alpha  level  is  accurate  asymptotically  in  the  number  of 
waves.  The  relevant  asymptotics  for  cross-lagged  panel  studies 
are  in  the  number  of  replications  since  the  number  of  waves  tends 
to  be  very  small. 

Mayer  (1984b)  considered  two  informal  procedures  for  detect¬ 
ing  the  presence  of  serial  correlation  in  the  errors  for  the 
multivariate  version  of  the  cross-lagged  panel  model.  His  first 
method  is  to  ignore  the  possibility  of  serial  correlation,  esti¬ 
mate  the  matrix  of  regression  coefficients  by  a  pooled  version  of 
least  squares  estimation,  and  then  use  the  residuals  to  estimate 
the  autoregressive  parameter  matrix  This  method  gives  a  crude 
estimate  of  £  but  does  not  produce  an  estimator  that  is  consistent 
as  the  number  of  replications  becomes  large.  His  second  method  is 
to  partition  the  matrix  ^  conformally  with  £  and  sec  the  off- 
diagonal  blocks  equal  to  £.  The  diagonal  blocks  of  £  are  estima¬ 
ted  by  expressing  each  subresponse  ve  tor  at  wave  t  as  a  function 
of  both  subresponse  vector,  at  wave  t-1.  The  other  subresponse  is 
considered  an  exogenous  predictor  and  a  multivariate  version  of  a 
standard  method  of  econometrics  is  used  to  estimate  the  diagonal 
block  of  j This  method  is  crude  in  that  ignoring  the  off- 
diagonal  blocks  of  £  may  seriously  bias  the  estimates  of  the  off- 
diagonal  elements  of  ]=!,  and  these  blocks  are  a  major  focus  of  Che 
regression  approach. 

3.  TESTING  FOR  SERIAL  CORRELATION 

We  begin  with  a  conditional  approach  to  testing  for  the 
presence  of  serial  correlation  in  the  error  structure  of  the 
multivariate  cross-lagged  panel  model,  conditional  in  the  sense 
that  it  treats  the  first  two  waves  as  fixed.  We  then  consider 
an  unconditional  approach  that  captures  the  random  nature  of  all 
the  waves.  We  complete  the  section  by  discussing  the  impacts  of 
serial  correlation  on  the  procedures  used  to  estimate  and  test 
the  matrix  of  regression  coefficients. 
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3 . 1  Conditional  Analysis 

Combining  equations  (2.1)  and  (2.2)  yields 

^t  "  ^c-1  $  +  ^c-1  ^  +  ^t  t“l,...,T 


Using  the  fact  chat  ^  -  ^_2  $ 

lt  "  3t-l(S  +  V  "  St-2  ($  *  It  C“2’--*’T  (3.1) 

which  expresses  the  cross-lagged  panel  model  with  first-order 
autoregressive  errors  as  a  second-order  vector-valued  autoregres¬ 
sive  model  with  independent  errors  over  waves  and  independent 
replications  at  each  wave,  a  model  studied  by  Anderson  (1978). 

We  suggest  chat  the  test  for  the  presence  of  serial 
correlation  be  preceded  by  a  test  for  no  effects  or  independence 
of  the  response  vector  over  time,  the  null  hypothesis,  ^ 

■  S£.  The  null  hypothesis  implies  H*:  £  +  £  -  £  and  ]JjJ  ■  jj,  a 
hypothesis  about  the  models  as  expressed  in  (3.1).  The  model  is 
reexpressed  as 


£  “  £-1  $  +  l 


where  £  -  ^ 


z *  •  .7'  r7ci)' 

^T-L  *  *  *  *  *xi  *~1 


-(2)' 

*-l 


£  *  '  and  $ 


We  define  the  estimators 


5  -  tCi  s^r1 


i +  £  K 


4  -  (n(T  -  digs  -  i_x  qrcs  -  ^  i 


and  Che  cesc  statistic 


C1  •  «$’  Ci  Cl  « k~l 

A  proof  quite  similar  to  those  in  Anderson  (1978)  yields 

Theorem  1:  Under  the  null  hypothesis  of  no  effects  (H*)  and 
given  and  ^ ,  the  test  statistic  is  a  Lawley-Hotelling  trace 
statistic  and  has  a  chi-square  distribution  vith  2k^  degrees  of 
freedom  as  the  number  of  replications  becomes  large. 

If  H*  is  rejected  two  possibilities  are  considered.  Either  £ 
and  ^  are  both  not  zero  or  £  alone  is  not  zero.  Because  of  the 
nature  of  the  model  we  dismiss  the  possibility  of  serially 
correlated  errors  (Q.  f  in  a  model  with  no  regression  effects 
(JJ  ■  JJ)  and  the  hypothesis  of  "countering  effects",  (JjJ  ■  -£) . 

We  differentiate  between  the  two  possibilities  by  testing  the 
hypothesis  H^:  ^  ■  £  which  is  equivalent  to  testing  the  hypothe¬ 
sis  H*:  A  ■  £  (since  we  concluded  £  +  £  as  a  result  of  the 

AAA  A 

first  test).  Let  ■  (QJ,  Qp  '  where  is  partitioned  conformally 
with 

ut  Ci  • 


■ «  k  ^ Vi’  <£y&-1  «£i IP’iy-1 


We  have 

Theorem  2:  Under  the  null  hypothesis  of  no  serial  correla¬ 
tion  (H*)  and  given  ^  and  the  test  statistic  t^  is  a  Lawley- 

Hotelling  trace  statistic  and  has  a  chi-square  distribution  with 
2 

k  degrees  of  freedom,  as  the  number  of  replications  becomes 
large. 

The  proof  of  this  result  is  a  direct  application  of  a  result 
of  Anderson  (1978). 

The  calculation  for  this  conditional  analysis,  can  be  done 
by  a  standard  statistical  package  such  as  SAS.  Two  limitations 


of  this  conditional  analysis  are  that  it  sacrifices  the  2nk 
degrees  of  freedom  of  the  first  two  waves  and  that  it  provides 
no  direct  estimate  of  the  regression  matrix  B. 

3.2  Unconditional  Analysis 

Applying  the  same  reasoning  as  used  in  the  conditional 
analysis,  we  begin  the  unconditional  analysis  by  considering  the 
hypothesis  of  no  effects  3  m  £  ■  0. 

The  likelihood  function  L(B,  A,  ©q,  generated  by  the 
nk(T  +  1)  observations  is 

nk(T  +  1)  _  nCT  -  1)  _n  _n 

«»>  2  1(1  2  If  lo4  *  4  I  2  l«ol  2 

h  £  ih  *  ‘r<f  -  *  f'1  -W* 


'  Sc-i($,  *  f  -  Zc-2(-W>f  'Sc  -  St-i'S  *  f  - 


(3.2) 


which  can  be  maximized  numerically. 

Under  the  null  hypothesis  the  likelihood  function 

L(0,  0,  A,  9- ,  0)  is 

•v.  \  x  \0  \ 

_  nk  (T  1)  _  nT  _  n 

«•>  2  If  2  l&l  2  «*p(-T!tr5o  &l 


t»l 

*  -1  *  -1 
which  is  maximized  by  letting  ^  ■  n  ^  “  (nT) 

T 

Z  V  2  . 
t-1  "C  "C 


(3.3) 


If  we  let  L  be  the  maximum  of  the  likelihood  in  (3.2)  and  L 
u  r 

be  the  maximum  of  the  likelihood  in  (3.3)  and  define 


L 

C.  »  -2  log  — 

L 

u 

we  have 

Theorem  3:  Under  Che  hypothesis  of  no  effect  ^ 

2 

the  test  statistic  t^  has  a  chi-square  distribution  with  (5k  ♦ 

k)/2  degrees  of  freedom  as  the  number  of  replications  becomes 
large. 

If  we  reject  ^  we  proceed  to  test  the  hypothesis 

of  no  serial  correlation  in  the  error  structure,  ^  by 

applying  the  likelihood  ratio  method.  Under  the  likelihood 
function  LQJ,  j\,  £)  is 

nk  (T  »  1)  _  nT  _  n 

(2«)  2  UJ  2  1^!  2  * 

«  r  <?t  -  <3-4> 
e-1 

Applying  a  result  of  Anderson  (1978)  the  likelihood  function 
is  maximized  by 

„  ,  T-l  T 

i  -  <nT>  t  2  Vi  2  *i-i  ?! 

i-0  i*L  11  1 

%  m  a  ^0  ^0 

and 

"i* 

t  «.!-?!_!  5)  '«.!-?!-!  ?> 

i»l 

Let  L*  be  the  maximum  of  the  likelihood  function  (3.4)  and  define 
r 

L* 

C  -  -2  log  — 

L 

u 

where  L  is  given  in  the  definition  of  t,:  we  have 
u  3 
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Theorem  4 :  Under  Che  hypothesis  of  no  serial  correlation 
(H^:  ^  ■  (?)  the  test  statistic  t^  has  a  chi-square  distribution 

with  (3k“  ♦  k)/2  degrees  of  freedom  as  the  number  of  replications 
becomes  large. 

The  tests  given  in  Theorem  3  and  Theorem  4  are  somewhat 
difficult  to  apply  in  that  both  require  maximization  of  the 
likelihood  function  given  in  (3.2),  a  calculation  which  requires 
,fbrute  force"  numerical  optimization.  With  a  little  modification, 
however,  tests  similar  to  those  in  Theorem  3  and  Theorem  4,  but 
more  easily  computed,  are  provided. 

Noce  that  the  maximization  of  the  likelihood  func'lon  given 
in  (3.2)  yields  maximum  likelihood  estimators  of  jj ,  and 

3  *  £  which  are  mathematically  consistent.  More  explicitly,  if 
we  let  JI,  ,  ,  and  be  these  estimators,  then  they  satisfy 

the  constraint 


(£U)  -  93  -  -3r  (3.5) 

If  we  relax  this  constraint  then  the  parameters  3,  3  +  1,  and 
34  can  be  estimated  independently.  This  relaxation  can  be 
achieved  by  replacing  the  model  of  the  wave  one  response  matrix  in 
(2.1)  by  the  model 


2  *2  3  -  r 

*4  %  *4 


(3.6) 


where  3q  Is  not  restricted  to  be  the  same  as  p. 

A  standard  maximum  likelihood  routine,  such  as  LISREL,  can 
be  used  to  maximize  the  likelihood  function  for  this  relaxed 
model.  Let  L(3,  3n ,  A,  A.  be  the  modified  likelihood 

->  A J  ,  -  -rf  J 

function  and  L*  m  L(3*.  3? »  A*,  A*,  3*,  J*  )its  maximum. 

U  w  '  *  Aj  v 

The  hypothesis  of  independence  of  the  response  vector  £ 
over  time  for  the  relaxed  model  is  ^  ^  “  S?  •  We  define 

the  test  statistic  for  H,  to  be 


L 

t*  *  -2  log  _ r  . 

J  L* 

u 
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For  the  relaxed  model  the  hypothesis  of  ao  serial  correlation 
remains  H^:  ^  ■  jj.  We  let  be  the  maximum  of  L(£,  ^ ,  jj),  j^, 
£)  and  define 

t*  ■  -2  log  **  r  . 

L* 

u 

By  the  usual  theory  of  likelihood  ratio  tests  we  have 

Theorem  5:  For  the  relaxed  model  as  the  number  of  replica¬ 
tions  becomes  large  under  5^  the  test  statistic  t*  has  a  chi- 

2 

square  distribution  with  (7k  +  k)/2  degrees  of  freedom,  and  under 

H.  the  test  statistic  t*  has  a  chi-square  distribution  with 

1  2  4 
(3k  +  k)/2  degrees  of  freedom. 

Note  that  under  the  relaxed  model  and  the  unconditional 

model  are  identical.  We  also  remark  that  although  the  hypothesis 

of  no  serial  correlation  is  the  same  for  the  two  models,  it  is  not 

obvious  how  the  difference  between  the  models  is  reflected  in 

hypothesis  tests. 


3.3  STATISTICAL  INFERENCES  WITH  SERIAL  CORRELATION 

Before  applying  our  tests  to  a  panel  model,  we  consider  the 
impact  of  serial  correlation  in  the  error  structure  on  the 
estimates  of  the  regression  parameters  and  the  tests  of  the 
hypotheses  of  no  effects  and  no  cross-effects. 

Mayer  (1984b)  has  shown  that  ignoring  serial  correlation 
leads  to  estimates  of  the  regression  parameters  that  do  not 
converge  to  the  true  parameters  as  the  number  of  replications 
becomes  large.  Similarly  the  likelihood  ratio  tests  for  the  test 
of  no  effects  and  no  cross-effects  can  be  highly  biased  if  serial 
correlation  is  ignored. 

Suppose  the  conditional  analysis  of  the  model  is  used  and 
serial  correlation  is  detected.  Then  the  parameter  matrices 
£  ♦  ^  are  estimable  but  the  regression  parameter  matrix  is 

not.  The  hypotheses  of  no  effects  and  no  cross-effects  can  be 
tested  by  testing  the  corresponding  hypothesis  about  ^  +  A  and  BA . 


If  Che  unconditional  analysis  is  used  and  che  constraints 
given  in  (3.5)  are  sacisfied  by  che  estimates  of  £  and  then 
expressions  are  given  for  £  and  £.  If  che  relaxed  model  is  used  - 
the  constraints  are  ignored  -  then  che  values  of  +  £  and 

can  be  used  to  solve  for  £  in  any  of  of  three  ways.  The  three 
solutions  are  only  asymptotically  identical. 

Whether  che  constraints  are  ignored  or  not  che  likelihood 
ratio  method  can  be  applied  to  give  tests  for  che  hypotheses  of  no 
effects  and  cross-effects. 

Several  other  methods  have  been  suggested  to  estimate  che 
regression  parameters  in  che  cross-lagged  panel  model  when  serial 
correlation  is  present  in  the  errors  (e.g.,  Markus  (1979),  Kessler 
and  Greenberg  (1980)).  One  method  is  to  ignore  the  off-diagonal 
elements  of  ^  and  to  express  each  response  at  time  t  as  a  linear 
function  of  che  responses  at  time  t-1  and  t-2  plus  an  error  term 
which  is  independent  over  time.  Then  ordinary  least  squares  is 
applied  to  that  equation  and  the  diagonal  elements  of  ^  are 
obtained  from  the  estimated  coefficients.  This  method  resembles  a 
standard  method  in  econometrics  attributed  to  Durbin  (Johnston 
(1972),  p.  263).  A  second  method  expresses  each  response  at  time 
t  as  a  function  of  che  responses  at  time  t-l  and  then  applies  che 
instrumental  variable  method  of  econometrics. 

Neither  of  these  methods  associates  naturally  with  a  test  of 
the  hypothesis  of  no  cross-effects.  Furthermore,  by  ignoring  the 
off-diagonal  elements  of  Che  estimates  of  the  off-diagonal 
elements  of  {!  may  be  highly  biased,  and  these  are  the  coefficients 
that  are  used  to  indicate  cross-effects.  Consequently,  a  test 
based  on  che  asymptotic  distribution  of  these  estimators  may 
seriously  over  estimate  che  degree  of  cross-effects. 

4.  ANALYSIS  OF  PANEL  DATA 

We  illustrate  our  methods  by  applying  them  to  a  sec  of  panel  data 
that  are  taken  from  a  study  of  che  relationship  between  patients 
opinion  of  Che  concept  of  Health  Maintenance  Organizations  (HMO) 
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and  their  perception  of  the  quality  of  care  they  are  receiving 
from  the  Health  Maintenance  Organization  in  which  they  are 
enrolled.  For  the  purpose  of  demonstrating  our  results,  the 
critical  issue  is  whether  a  cross-lagged  panel  model  fit  to  these 
data  appears  to  have  serial  correlation  in  the  error  structure. 
Secondary  issues  include  the  degree  to  which  such  correlation 
affects  the  estimates  of  the  regression  parameters  and  the  degree 
to  which  the  simpler  conditional  analysis  produces  results  that 
are  similar  to  the  results  produced  by  the  unconditional  analysis. 

These  data  were  obtained  form  interviews  of  20  patients  of  a 
particular  HMO.  They  were  selected  at  random  and  Interviewed 
after  having  seen  a  primary-care  physician  at  the  HMO  for  the 
first  time.  The  same  20  patients  are  interviewed  again,  after 
their  second  and  third  visits  to  the  HMO.  The  two  variables  of 
interest  are  each  a  compilation  of  several  measurements  that  are 
made  at  each  interview.  The  original  interviews  were  not  provided 
by  the  consulting  firm  that  owns  the  data. 

The  first  variable  measures  the  degree  to  which  the  patient 
supports  the  concept  of  an  HMO  as  a  provider  of  primary  medical 
care.  Issues  of  secondary  and  emergency  care  were  not  addressed 
in  the  study.  The  second  variable  measures  the  degree  to  which 
the  patient  is  satisfied  with  the  treatment  received  at  the  HMO  on 
the  particular  visit  just  concluded.  The  variables  are  standard¬ 
ized  (over  a  larger  sample)  to  have  means  of  10  and  standard 
deviations  of  3.  The  panel  model  adopted  is  the  simplest  example 
of  the  multivariate  model  since  p  ■  q  ■  1. 

The  conditional  analysis  was  completed  by  the  multivariate 
regression  routine  in  SAS  and  produced  the  results  in  Table  1. 
Displayed  are  the  the  Lawley-Ho celling  trace  statistic  for  testing 
the  hypothesis  of  no  effects,  c^  ■  342,  which  is  asymptotically 
chi-square  with  8  degrees  of  freedom  and  the  Lawley-Ho telling 
trace  statistic  for  testing  the  hypothesis  of  no  serial  correla¬ 
tion  in  the  error  structure,  t^  *  31.1,  which  is  asymptotically 
chi-square  with  4  degrees  of  freedom.  Combining  these  tests  gives 
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TABLE  1 


CONDITIONAL  ANALYSIS  OF  HMO  PANEL  DATA 


Maximum  Likelihood  Estimates 


= 


l*! 


* : 


,715 

.041' 

estimates 

*  +  k 

,149 

.591 

,338  ■ 

-.137 

*  H 

,645 

.548 

|  estimates 

i  Likelihood 

Estimates  of 

Standard  Errors 

,080 

.116  ^ 

,061 

.088 

,089 

.1291 

,154 

.224 

of  the 

Estimates  to  Their  Estimated  Standard 

9.0 

.35] 

-2.4 

6'7J 

3.8 

-1.06  ! 

-4.2 

2.45 

Lavley-Ho telling  Test  of  H?:  £  -  £  and  ££  *  £ 

asymptotic  chi-square  8df 
Lawley-Hotelling  Test  of  H|:  ^  ■  £ 
t^  »  31. L  asymptotic  chi-square  4df 


tL  -  342 
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evidence  Chat  there  are  temporal  effects  both  in  the  observed 
regression  and  in  the  error  structure. 

Also  displayed  in  Table  1  are  the  maximum  likelihood 
estimates  of  the  elements  in  £  +  £  and  the  corresponding 
estimates  of  the  standard  errors  of  these  elements,  and  the  ratios 
of  the  estimated  coefficients  to  their  estimated  standard  errors. 


Examination  of  these  coefficients  and  the  standard  errors 


indicates  responses  at  a  giver,  wave  are  good  predictors  of  the 
same  responses  at  the  next  wave.  Also  indicated  is  that  the 
patient's  perception  of  the  quality  of  care  at  wave  t  may  have  a 
stronger  affect  on  his  or  her  opinion  of  the  concept  of  an  HMD  at 
wave  t  +  1  than  does  his  opinion  of  the  concept  of  an  HMO  at  wave 
t  have  on  his  perception  of  the  quality  of  his  care  at  wave  t  +  1. 

Note  again  that  this  conditional  analysis  provides  no  direct 
estimate  of  the  regression  parameter  matrix  or  autoregressive 


parameter  matrix. 

Before  turning  to  the  unconditional  analysis  we  examine  Table 


2  which  displays  the  results  that  would  have  been  obtained  if  the 
conditional  analysis  were  applied  to  the  model  without  allowing 
for  serial  correlation  in  the  errors.  Remembering  that  the  model 
being  analyzed  is  incorrectly  specified,  we  note  that,  the 
Lawley-Ho telling  trace  statistic  for  testing  the  significance  of 
the  regression  parameters  (22.9),  which  is  "advertised"  to  be 
asymptotically  chi-square  with  4  degrees  of  freedom,  is  quite  a 


bit  smaller  than  the  test  statistic  t^  given  in  Table  1. 


Examination  of  the  estimated  coefficients  and  the  estimated 


standard  errors  Indicates  far  less  temporal  dependence  than  is 
indicated  by  the  results  in  Table  1.  Furthermore,  these  results 
appear  to  indicate  that  the  patients  opinion  of  the  concept  of  an 
HMO  at  wave  t  may  have  a  stronger  impact  on  his  or  her  perception 
of  the  quality  of  care  received  at  wave  t  +  1  than  does  his  or  her 
perception  of  the  quality  of  care  received  at  wave  t  has  on  his  or 
her  opinion  of  the  concept  of  HMO  at  wave  t  +  1;  this  result  is  in 
sharp  contrast  to  the  results  of  Table  1  and  indicates  that 


TABLE  2 

CONDITIONAL  ANALYSIS  OF  HMO  PANEL  DATA: 
IGNORING  SERIAL  CORRELATION 


Maximum  Likelihood  Estimates 


529 

-.313 

123 

.360 

Maximum  Likelihood  Estimates  of  Standard  Errors 
.153" 

.155 

Estimates  to  Their  Estimated  Standard  Errors 

-2.o" 

2.3 

"Lawley-Hotelling"  test  of  H^.  ^  [Assuming  £  »  £] 
t  ■  22.9  asymptotic  chi-square  4df 

ignoring  the  serial  correlation  in  error  structure  would  have  led 
to  misleading  conclusions. 

The  unconditional  analysis  of  the  relaxed  model  was  completed 
with  the  LISREL  program  associated  with  SPSS-X  and  produced  the 
statistics  presented  in  Table  3.  This  table  gives  the  approximate 
likelihood  ratio  statistic  for  testing  the  fit  of  the  model  when 
compared  to  the  model  of  no  effects,  t*  -  86.81,  which  is 
asymptotically  chi-square  with  15  degrees  of  freedom  and  the 
approximate  likelihood  ratio  statistics  for  testing  for  the 
presence  of  serial  correlation  in  the  error  structure,  t*  - 
47.76,  which  is  asymptotically  chi-square  with  7  degrees  of 
freedom.  The  results  of  the  tests  are  consistent  with  the  results 


l  : 


.138 

.140 


Ratios  of  the 
3.8 
I -.88 


* ! 


A'.V.N 


TABLE  3 


UNCONDITIONAL  ANALYSIS  OF  HMO  PANEL  DATA: 
RELAXED  MODEL 
(Continued) 


’  9.54 

.38' 

It*' 

2.91 

-.871 

-2.61 

7.11 

-.871 

2.91 

Likelihood  Ratio  Test  of  ^ 

t*  *  86.81  asymptotic  chi-square  15df 

Likelihood  Ratio  Test  of  H^:  ^  -  Q 
t*  ■  47.76  asymptotic  chi-square  7df 

of  the  tests  from  the  conditional  analysis.  We  conclude  that  the 
panel  model  has  both  regression  effects  and  serial  correlation  in 
the  error  structure. 

Also  displayed  in  Table  3  are  the  maximum  likelihood  esti¬ 
mates  of  the  parameters,  their  estimated  standard  errors  and  the 
ratio  of  the  estimates  to  the  estimated  standard  error.  These 
statistics  support  conclusions  almost  identical  to  those  obtained 
from  examination  of  the  statistics  obtained  from  the  conditional 
analysis  (Table  1).  The  estimates  of  ^  ♦  £  and  are  similar 
although  the  estimates  of  the  standard  errors  differ.  We  suggest 
that  given  the  ease  of  calculation  and  the  similarity  of  results 
the  conditional  analysis  might  be  more  attractive  to  the  practi¬ 
tioner.  Table  4  displays  the  unconditional  analysis  of  the  model 
with  the  serial  correlation  Ignored.  We  note  that  these  results, 
like  the  results  in  Table  2,  could  be  quite  misleading.  Ignoring 
the  serial  correlation  in  the  error  of  a  panel  model  appears  to  be 
a  serious  error  whether  the  model  is  treated  conditionally  or 
unconditionally. 
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TABLE  4 


UNCONDITIONAL  ANALYSIS  OF  HMO  PANEL  DATA: 
IGNORING  SERIAL  CORRELATION 


(Nominal)  Maximum  Likelihood  Escimaces 


r 


t* 


.536 

-.298 

-.106 

.401 

’  5.57 

-1.90 

-1.90 

6.66 

*7.657 

2.911 

2.911 

2.983 

(Nominal)  Maximum  Likelihood  Estimates  of  Standard  Errors 

r* 

Ratio  of  the  Estimates  to  Their  Estimated  Standard  Errors 


.147 

.161 

.180 

.197 

1.350 

1.093 

1.093 

1.615 

r 


t*s 


3.6  -1.8 

-.589  2.036 

'4.123  -1.734 

-1.734  4.123 


Likelihood  Test  of  Hq:  £  ■  £  [Assuming  £  »  £] 
19.72  asymptotic  chi-square  4df 


5. _ DISCUSSION 


We  have  developed  and  demonstrated  tests  for  detecting  the 
presence  of  serial  correlation  in  the  error  structure  of  a  cross- 
lagged  panel  study.  Allowing  the  error  structure  to  contain 
serial  correlation  may  extend  the  usefulness  of  the  multivariate 
regression  model  since  observations  made  on  a  single  subject  at 
two  waves  are  rarely  independent  particularly  if  the  waves  are 
close. 

A  second  method  of  allowing  dependence  in  the  error  structure 
has  been  formulated  by  econometricians  for  the  univariate  continu¬ 
ous  variable  panel  model  (eg,  Balestra  and  Nerlove  (1966),  and  Wallace 
and  Hussain  (1969)).  In  this  formulation  the  error  for  a  given 
observation  can  be  decomposed  into  the  sum  of  a  pure  error  and 
replication  effect.  The  replication  effect  captures  the  tendency  for 
sampling  units  (subjects)  that  are  above  the  regression  line  at  the 
first  wave  to  stay  above  the  regression  line  across  waves.  Our  work 
with  business  data  has  convinced  us  that  this  error  structure  is  often 
more  realistic  than  the  independent  error  structure  and  is  a  viable 
competitor  to  the  serial  correlated  structure. 

Anderson  and  Tsiao  (1981;  1982)  have  studied  the  problem  of 
estimating  the  parameters  of  the  univariate  panel  model  with  this 
decomposable  error.  We  are  currently  extending  their  results  to 
the  multivariate  model  and  considering  the  problem  of  distinguish¬ 
ing  between  errors  that  are  serially  correlated  and  errors  that 
are  decomposable. 


*The  authors  gratefully  acknowledge  discussions  with 
colleagues  and  students.  Particular  thanks  go  to  T.  W.  Anderson 
and  D.  R.  Rogosa. 
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